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ABSTRACT 

We use state-of-the-art radiation-MHD simulations and 3D non-LTE radiative transfer computations 
to investigate Ha line formation in the solar chromosphere and apply the results of this investigation 
to develop the potential of Ha as diagnostic of the chromosphere. 

We show that one can accurately model Ha line formation assuming statistical equilibrium and 
complete frequency redistribution provided the computation of the model atmosphere included non- 
equilibrium ionization of hydrogen, and the Lyman-a and Lyman-/3 line profiles are described by 
Doppler profiles. 

We find that 3D radiative transfer is essential in modeling hydrogen lines due to the low photon 
destruction probability in Ha. The Ha opacity in the upper chromosphere is mainly sensitive to the 
mass density and only weakly sensitive to temperature. 

We find that the Ha line-core intensity is correlated with the average formation height: the larger 
the average formation height, the lower the intensity. The line-core width is a measure of the gas 
temperature in the line-forming region. The fibril-like dark structures seen in Ha line-core images 
computed from our model atmosphere are tracing magnetic field lines. These structures are caused by 
field-aligned ridges of enhanced chromospheric mass density that raise their average formation height, 
and therefore makes them appear dark against their deeper- formed surroundings. We compare with 
observations, and find that the simulated line-core widths are very similar to the observed ones, 
without the need for additional microturbulence. 

Subject headings: Sun: atmosphere — Sun: chromosphere — radiative transfer — magnetohydrody- 
namics (MHD) 



1. INTRODUCTION 

The Ha line is one of the most popular lines for study- 
ing the solar chromosphere. In fact, the chromosphere 
was defined as what is seen in this line (Lockver 186^. 
Yet, nearly 150 years after its discovery there is still much 
unknown about the chromosphere and the formation of 
Ha. For an overview of the current state of understand- 
ing and open questions of Ha as a diagnostic of the chro- 
mosphere we r efer to a n exce llent series of discussions by 
IRuttenl (l200l,|20T0a,b, 2011). 

The lack of understanding is due to the particularly 
complicated physics of this layer. It forms the inter- 
face between the convection zone and the corona and 
chromospheric modeling therefore has to include at least 
the upper convection zone and the lower corona. The 
presence of magnetic fields leads to structures that can- 
not be reproduced in one or two-dimensional modeling. 
Three-dimensional geometry is therefore required, which 
significantly increases the amount of computational work 
required. 

Over the last few years it has become possible to cre- 
ate models of the chromosphere that include much of 
the required physics and the required physical extent: 
they model the atmosphere from the upper convection 
zone to the lower corona, and include resistive magne- 
tohydrodynamics, full 3D radiation in the photosphere 
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and lower chromosphere, parameterized radiative losses 
in the upper chromosphere and corona, the rmal con- 
duction and a realistic equation of state. (iHansteenl 
I200I iHalsteen et al.l 120071: IGudiksen et al.ll2011[ ) These 
models have been used to study various chromos pheric 
phenomena, such as dynam ic fibrils (Hansteen et al.l 
120061: iDe Pontieu et al.ll2007aD. t he presence of Alfvenic 
waves jibe Pontieu et al.l l2007b D_, r nagnet ic flux emer- 
gence (. Martfn ez-Svkora e t al. 2 0081 [2009al ) and spicules 
(IMartfnez-Svkora et al-llkOObt i20ll 

All these studies compared the time-evolution of the 
MHD variables in the simulations with some observed 
quantity, rather than first computing synthetic observa- 
tions. For example, the dynamic fibril studies compared 
observed motions of structures observed in Ha to mo- 
tions of a temperature iso-contour in the numerical sim- 
ulations. 

Direct comparison of synthetic observables generated 
from simulations with observations have been more 
scarce. This is partly caused by the complex radia- 
tive transfer processes in the chromosphere. Spectral 
lines form in non-LTE, partial redistribution of photons 
(PRD) is important for many lines, and the effects of 
full 3D radiation require complicated and computing- 
intensive radiative transfer calculations. Examples of 
such computations are the work on bisectors of the Na I D 
and Ca II 8542 A lines bv lUitenbroeg ()2006bf ) , who used a 
3D hydrodynamical model that extended up to 1000 km 
above optical depth unity, covering the lower and mid 
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chromosphere only. Leenaarts fc Carlssonl ()2009[ ) and 
iLeenaarts et"all ()200a |2010[ )' investigated the formation 
of Na ID, Ca II 8542 A and Ca II H in 3D MHD models 
that extended up into the corona. They found that 3D 
effects are important in the cores of all lines. The 3D 
effects reduced the overall contrast of line core images, 
and for Na I D it was found that the bright halos around 
magnetic field concentrations are caused by photon scat- 
tering. 

The previous studies found reasonable agreement be- 
tween observations and synthetic imagery for the low- 
forming Na ID line. Agreement for the Ca II 8542 A was 
reasonable for qui et reg ions without overlying fibrils. 

As discussed bv ILeen aarts (20 i3) the discrepancy be- 
tween observed and synthetic Ha was big: computa- 
tions treating each column as a plane-parallel atmosphere 
(henceforth ID radiative transfer) showed a granulation 
pattern in the line core. This is in stark contrast to 
the observations th at show either fibrils or the signature 
of acoustic shocks (jRutten et al.l 120081 ). It was specu- 
lated that this discrepancy could be caused by a lack 
of spatial extent and spatial resolution of the employed 
simulations, lack of treatment of non-equilibrium level 
populations (with non-equilibrium it is meant that the 
time derivative of the level populations is not necessarily 
equal to zero) effects and/or the neglect of PRD. 

In the current paper we report an investigation of the 
formation of the Ha line. In Sec. [2HS] we show that in 
order to accurately model Ha line formation one has to 
use a model atmosphere computed with an equation of 
state (EOS) that takes into account the effect of non- 
equilibrium hydrogen ionization and one has to compute 
the radiation field in full 3D. We also show that one can 
neglect the time-dependence of the rate equations and 
use statistical equilibrium in the subsequent non-LTE 
spectrum synthesis and that one can approximate PRD 
effects by treating the Lyman lines with Doppler profiles 
in complete redistribution. In Sec. |6H8] we analyze the 
results of our 3D computation in detail, compare with 
observations and finish with a discussion of the results. 

2. THE EFFECT OF NON-EQUILIBRIUM TRANSFER ON 

Ha 

To investigate the effect of non-equilibrium radiative 
transfer on the formation of the Ha line we use a time- 
series of hydrodynamical snapshots computed with the 
RADYNcode (Carlsson & Stein 1992, 1995, 1997, 2002). 

This code solves the equations of conservation of mass, 
momentum, energy and mass. In addition it solves 
the non-LTE radiative transfer problem including non- 
equilibrium rate equations for hydrogen, helium and cal- 
cium. It uses a 5-level plus continuum model atom for 
H I and Ca II, and a nine-level helium atom with 5, 3 
and 1 levels in He I, He II and He III respectively. These 
equations are solved implicitly on an adaptive mesh that 
moves po ints to regions where high spatial resolution is 
required (Dorfi fc Drurvlll987D . RADYN solves the non- 
equilibrium radiative transfer problem for hydrogen self- 
consistently including its effect on the ionization balance 
and therefore on the EOS. It treats PRD effects in the 
Lyman lines rather crudely, by truncating the Lyman 
line profiles at 48kms^^ away from the line-center fre- 
quency. We show in Sec. |4] that this simple method has 
only a minor effect on the formation of Ha. 




0.0 



l).5 1.0 
z [Mm| 



1.5 



2.(1 




Figure 1. Comparison of non-equilibrium radiative transfer (red 
curves) to radiative transfer assuming statistical equilibrium (blue 
curves) in a one-dimensional model atmosphere computed with Ra- 
dyn. Top: level population densities of the lower (solid) and upper 
(dashed) level of the Ha line as function of height. Middle: gas 
temperature (black) and the Ho line source function as function 
of height. The additional blue and red scales show log^g r at line 
center. Bottom: Vertically emergent Ha line profile. 

In our RADYN simulation the lower boundary is lo- 
cated just below the photosphere at a column mass of 
8.1gcm~^, roughly corresponding to rsoonm = 16. The 
upper boundary is located in the corona, at 10,000 km 
above the lower boundary. The bottom boundary is 
closed and driven with a prescribed velocity field derived 
from MDI observations of the Doppler shifts of the 676.78 
nm Ni I line in a quiet sun region. This driver excites 
waves that travel upward and steepen into shocks. We 
ran the simulation for slightly more than 1 hour of solar 
time and saved output at 10 s intervals, yielding a series 
of 374 snapshots of both the hydrodynamical state and 
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Figure 2. Comparison of non-equilibrium radiative transfer (red 
curves) to radiative transfer assuming statistical equilibrium (blue 
curves) in a time-series of one-dimensional model atmospheres com- 
puted with Radyn. Top: vertically emergent intensity in the Ha 
line core. Bottom: average formation height. 



the pertinent radiative transfer quantities. 

For each of these ID snapshots we then solved the non- 
LTE radiative transfer problem for hydrogen assuming 
statistica l equilibr i um w ith the radiative transfer code 
MULTI (jCarlssonI Il986n using identical hydrodynamic 
quantities (temperature, mass density, electron density 
and vertical velocity) and model atom. MULTI and RA- 
DYN have identical source code for the routines that 
compute the radiative transfer. We thus ensured that 
any differences in the results are caused by the change 
from non-equilibrium (NE) level populations to transfer 
assuming statistical equilibrium (SE), and not by differ- 
ences in the numerical method. 

Figure [T] compares the formation of the Ha line be- 
tween NE and SE calculations a,t t = 3740 s. At this 
time there is a strong shock in the chromosphere, a cir- 
cumstance which typically exhibits the largest differences 
in line formation between NE and SE. At times that the 
chromospheric temperature structure is smooth there is 
almost no difference between the two assumptions. 

The top panel shows the population densities of its up- 
per and lower level. In SE both populations are slightly 
higher than in the NE case. Their ratio is however 
nearly equal, leading to a very similar source function, as 
shown in the middle panel. The SE computation shows 
a slightly higher source function throughout the simu- 
lated chromosphere. The line-core optical depth scales 
are nearly equal. The bottom panel compares the line 
profiles. They are near-identical, except in the line core, 
where the SE intensity is larger, which reflects the higher 
source function. 

In Fig. [2] we compare the time evolution of the line- 
center intensity and the line-center average formation 
height Za- The latter is defined as 



(1) 



where Zh and Zt are the height of the bottom and top 
boundary in the model and x, S and r are the line-center 
opacity, source function and optical depth, respectively. 
The average formation height is in general not a good 
measure of where the Ha line-core intensity is formed; 
the contribution functions are often multiply peaked and 
can be non-zero over ranges of more than 1 Mm. How- 
ever, if the formation height range has to be quantified 
by one number, the average formation height is a better 
measure than the height of optical depth unity. 

The curves for SE closely follow the corresponding 
curves for the full non-equilibrium solution. The SE 
intensity is consistently higher than the NE intensity, 
varying between 10% higher in the presence of shocks 
to only 2% higher during inter-shock phases. The SE 
average formation height is typically 30 km higher dur- 
ing inter-shock phases, and 80 km higher during shock 
passages. Both the differences in intensity and average 
formation height are much smaller than the intrinsic vari- 
ability caused by the changing state of the atmosphere. 

Therefore, we conclude that it is possible to assume SE 
for the radiative transfer of hydrogen, provided the tem- 
perature and electron density in the model atmosphere 
are computed with a simulation that includes the effect 
of non-equilibrium hydrogen ionization. The latter is 
essential: assuming statistical equilibrium when calcu- 
lating the atmosphere leads to temperatures that can 
be more than 2000 K off and ele ctron densities that can 
be even orders of magnitude off (iCarlsson fc SteinI 120021 : 
iLeenaarts et al., . 20071 . 

This result is crucial in making it possible to compute 
realistic synthetic imagery in hydrogen lines from 3D 
simulations of the solar chromosphere. The detailed non- 
equilibrium treatment of radiative transfer as in RADYN 
is too computationally demanding in multi-dimensional 
geometry. It is, however, possible to compute the non- 
equilibrium ionization balance of hydrogen and its effect 
on the equation-of-state using approximations that are 
efficient enough to be included in a 3D MHD code. 

3. 3D SIMULATION SNAPSHOT AND 3D RADIATIVE 
TRANSFER 

In Sec. m and [5] we investigate the influence of partial 
redistribution in Lyman lines and 3D evaluation of the 
radiation field on the formation of Ha. We study these 
effects in a snapshot of a 3D radiation-MHD simulation 
computed with the Bifrost code. 

Bifrost solves the equations of resistive MHD on a stag- 
gered Cartesian grid. In addition to this the code can 
solve various additional physical processes. The snapshot 
we are using was computed with the following additional 
physics: Optically thick radiative transfer in the photo- 
sphere and low chromosphere, parameterized radiative 
losses in the upper chromosphere, transition region and 
corona, thermal conduction along magnetic field lines 
and EOS that includes the effects of non-equilibrium ion- 
ization of hydrogen. As discussed in Sec. [2l the latter is 
of critical importance for the subsequent computation of 
solution of the radiative transfer problem for hydrogen. 

This non-equilibrium EOS essentially solves the non- 
equilibrium rate equations for hydrogen with prescribed 
radiative rates, together with equations for energy, par- 
ticle and charge conservation. By using prescribed rates 
there is no need to solve the transfer equation, and there- 
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Figure 3. Characterization of the Bifrost snapshot. Upper left: vertically emergent intensity in the continuum radiation bin. Upper right: 
Vertical magnetic field at z = OMm. Lower left: Gas temperature at z = 1.7Mm. Lower right: Angle of the magnetic field vector with the 
horizontal, with the red arrows indicating the direction of the horizontal component of the magnetic field. 



fore the problem becomes computationally tractable in 
3D simulations. 

We refer the reader to iLeenaarts et al.l (|2007[ ) for 
more details about t he non-equilibrium EOS, and to 
IGudiksen et al.l ()2011D . for an extensive description of 
Bifrost. 

Our simulation has a grid of 504 x 504 x 496 grid 
cells, with an extent of 24 x 24 x 16.8 Mm. In the ver- 
tical direction the grid extends from 2.4 Mm below to 
14.4 Mm above average optical depth unity at 500 nm, 
encompassing the upper convection zone, photosphere, 
chromosphere and lower corona. The x and y-axes are 
equidistant with a grid spacing of 48 km. The z-axis is 
non-equidistant. It has a grid spacing of 19 km between 
z — —I and z — 5 Mm, while the spacing increases to- 
wards the upper and lower boundaries to a maximum 
of 98 km at the coronal boundary. The magnetic field 
has a predominantly bipolar structure, that manifests 
itself in the photosphere as two clusters of magnetic 
concentrations of opposite polarity. The magnetic field 
was introduced into a relaxed hydrodynamical simulation 
through specifying the vertical magnetic field at the bot- 
tom boundary with a potential field extrapolation giving 



the magnetic field throughout the box. The field at the 
bottom boundary was given with an averaged signed field 
of zero and two patches of opposite polarity separated by 
8 Mm. The simulation ran for 3000 s before switching on 
hydrogen non-equilibrium ionization. The snapshot we 
analyze is after 850 s with hydrogen non-equilibrium ion- 
ization. 

The average unsigned magnetic field strength at z = 
OMm in the snapshot we selected is 50 G. The snapshot 
is further characterized in Fig. [S] The upper-left panel 
displays the vertically emergent intensity in the contin- 
uum radiation bin. It shows granules and some bright 
intergranular structures that correspond to kG magnetic 
field concentrations as can be seen by comparing with 
the upper-right panel, which shows the vertical magnetic 
field strength a.t z — OMm. The two clusters of oppo- 
site polarity field stand out in black and white. The 
lower left panel shows the gas temperature in the chro- 
mosphere at a height of 1.7 Mm. The corona dips down 
to below this height above the photospheric magnetic 
field concentrations. Between the magnetic patches the 
temperature shows filamentary structure that follows the 
general orientation of the magnetic field. This can be 
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seen in the lower-right panel, that shows the angle of 
the magnetic field vector with the horizontal, together 
with arrows pointing in the direction of the horizontal 
component of the field, also at z = 1.7 Mm. 

We performed the radiative transfer computations 
based on this snapshot with the MultiSd code 
(|Leenaarts &: CarlssonI [20091 ) . This is a radiative trans- 
fer code capable of handling 3D model atmospheres. It 
can evaluate the radiation field for each column of the 
atmosphere separately (column-by-column plane-parallel 
approximation), or it can do so in full 3D, taking the hor- 
izontal structure into account. 

MultiSd employs t he accelerated lambda ite ration 
method developed by iRvbicki fc Humm"eil (|1992[ ) with 
the extension to treat effects of partial frequency re- 
distribution using the angle-averaged approximation by 
lUitenbroeld ()2001l ). In order to make the 3D radiative 
transfer problem computationally tractable we halved 
the horizontal resolution of the MHD snapshot and re- 
sampled the vertical grid to 200 equidistant points be- 
tween z — —0.5 Mm and z = 5.5 Mm to cover just the 
formation range of the various hydrogen transitions. For 
the radiative transfer computations we thus had a hor- 
izontal resolution of 96 km and a vertical resolution of 
32 km. 

We used a 5-level plus continuum hydrogen model 
atom. For computations assuming ID plane-parallel ge- 
ometry we used a 5-point Gauss-Legendre quadrature to 
compute the angle integration for each ray direction (so 
10 rays in total) for the radiation fie ld. For full 31) com- 
putations we employed the A4 set of lCarlsonl ()1963f ) with 
3 rays per octant, and thus 24 rays in total. MultiSd does 
not include an equation for charge conservation. The 
non-LTE proton density we obtained is larger than the 
electron density at certain locations in the model atmo- 
sphere. 

4. INFLUENCE OF PRD IN LYMAN LINES ON Ha 

The Lyman lines of hydrogen are strongly affected 
by partial redis tribution of photons over the line pro- 
file ( PRD, e.g.. iMilkev fc Mihalas|[l973t iVernazza et al.l 
|1981|) . The Lyman line transfer affects the atomic 
level populations and therefore also the formation of the 
Ha line. The computation of PRD effects in a three- 
dimensional atmosphere is currently unfeasible. This is 
partly due to the increase in computational cost com- 
pared to using the simplifying assumption of complete 
redistribution of photons over the line profile (CRD), 
and partly due to algorithm stability: the strong gra- 
dients present in MHD models very often lead to non- 
convergence of the A-iterations. We therefore investi- 
gated the effect of various approximations to PRD ef- 
fects in the Lyman lines, in the hope of finding a simple 
recipe that has the same influence on the formation Ha 
as PRD. 

Inspired by iMilkey fc MihalasI (|1973[ ) and 
lAvrett fc Vernazzal ( 19771 ) we compared four differ- 
ent treatments of the Lyman-a and /3 lines: using 
angle-averaged PRD with a Voigt profile (hereafter 
'PRD case'), CRD with a Voigt profile ('CRD case'), 
CRD with a Voigt profile truncated at 48 kms~^ away 
from line center ('truncated case') and CRD with a 
Gaussian profile with Doppler broadening only ('Doppler 
case'). 
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Figure 4. Comparison of Ha line formation between various 
methods to treat radiative transfer in the Lyman-o and /9 lines 
in a column of a 3D model atmosphere. In all panels the red 
curve indicate results obtained with angle-averaged PRD, orange 
with CRD, green with a Doppler profile and blue with a truncated 
Voigt profile. Top: optical depth versus height, the thin horizon- 
tal lines indicate r equals 0.1,1 and 10, to guide the eye. Middle: 
gas temperature (black curve) versus height, with the total source 
function overplotted. Bottom: contribution function to intensity of 
the Ha line center versus height. The thin horizontal line indicates 
dZ/dz = to guide the eye. The average formation height is indi- 
cated by the triangles. Emergent line profiles for this atmosphere 
are shown in Fig. \S\ 

We took an xz-slice from the 3D Bifrost snapshot at 
y = Mm and solved the SE radiative transfer problem 
for hydrogen with MultiSd for each of the four differ- 
ent approximations, but otherwise identical input data. 
Each of the 252 columns was treated as a plane-parallel 
atmosphere. It turned out that the differences between 
the various cases showed behavior common to nearly all 
columns. We therefore show the results for one represen- 
tative column only in Figs. |3]and[SJ 

The top panel of Fig. S] shows the optical depth as 
function of height for the various cases. The CRD case 
shows the lowest optical depth of the Ha core at any 
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Figure 5. Comparison of the Ha line profile computed with the 
Lyman-a and /3 lines treated with angle-averaged PRD (red), CRD 
(orange), a Doppler profile (green) and a truncated Voigt profile 
(blue). Details on how the line core is formed are shown in Fig.|4] 
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Figure 6. Comparison of the H« line core formation with the 
Lyman-a and /3 lines treated with angle-averaged PRD (red), CRD 
(orange), a Doppler profile (green) and a truncated Voigt profile 
(blue) in a 2D slice of a 3D radiation-MHD snapshot computed 
with Bifrost. Top panel: line core intensity versus spatial location. 
Bottom panel: average formation height versus spatial location. 
The stars indicate the column used in Figs. |4] and [5] 

given height, the truncated case the highest; the Doppler 
case is closest to the PRD case. In all cases there is no 
increase in optical depth with geometrical depth around 
z = 0.8 Mm, indicating the 'opacity gap ' typical of lines 
with an excited state as a lower level (jLeenaarts et al.l 
[20061) . This can be seen in observations when stepping 
through a series of Ha images at successively larger wave- 
length offset from line center, where the overall appear- 
ance jumps from high-chromospheric structure to photo- 
spheric granulation, without intermediate reversed gran- 
ulation such as in the Ca II H line. 

The ordering of the optical depth scales can be ex- 
plained in terms of the thermalization depths of the 
Lyman-a line, or in other terms, as the probability in 



Table 1 

Comparison of PRD approximations 



Case 


< I/IpRD - 1 > 


< - 2a, PRD > (km) 


Voigt 


0.026 


-203 


Truncated 


0.019 


96 


Doppler 


0.011 


65 



photon escape in the line wing fcf. lMihalaslfTOTl . The 
Doppler and truncated case have a very weak or no line 
wing, and therefore a low thermalization depth and few 
line photons escape through the line wings. In contrast, 
the Voigt case has prominent wings, a large thermaliza- 
tion depth and a fraction of the line photons escape into 
space through the wings. Therefore the Lyman-a line 
has a larger net downward rate and hence lower n — 2 
population and hence lower Ha opacity in the Voigt case 
than in the Doppler and truncated cases. The PRD case 
is intermediate: Lyman-a and Lyman-/? behave more like 
Doppler lines in the upper chromosphere and more like 
a Voigt lines deeper down in the atmosphere. 

The middle panel of Fig. 2] shows a comparison of the 
gas temperature and the total source function S in the 
Ha line core. The temperature drops down from the pho- 
tosphere to a minimum at z = 1.25 Mm, and then slowly 
increases to transition-region temperatures, with a dip 
at z = 2.0 Mm. All source functions show S = B (with 
B the Planck function) up to z = 0.2 Mm after which 
they decouple from the local temperature. The dips in 
the source functions at around z — 0.8 Mm are caused by 
the very low line opacity causing the background source 
function to dominate the total source function. The large 
spike at z = 2.5 Mm is likewise caused by the background 
processes and plays no role in the formation of the line. 
Except where the background source function is appre- 
ciable, in the core-forming region (0.3 < z < 2.5 Mm) the 
truncated and Doppler cases show a larger source func- 
tion than in PRD. The CRD case a lower source function. 
The bottom panel shows the contribution function (CF) 
to intensity. All CFs are multiple-peaked and very wide, 
being non-zero over a height range of 2 Mm. The CRD 
case has the lowest average formation height (za, indi- 
cated by a triangle), the truncated case the highest. 

Fig. [5] shows the vertically emergent line profile for the 
same column as shown in FigHl The profiles are identical 
in the line wing. In the ±0.6 A around line center the 
differences are small but visible. 

The top panel of Fig. [6] compares the variation of the 
line-core intensity and average formation height for all 
columns in the xz-slice. The variations of the Ha in- 
tensity due to the various PRD approximations in Ly-a 
and Ly-/3 are much smaller than the variations due to 
the different atmospheric properties. Table [T] compares 
the relative difference in line core intensity between the 
various approximations and the PRD case. The Doppler 
case gives the best approximation to the PRD case, with 
on average only 1.1% difference in line-core intensity. 

The bottom panel similarly compares the average for- 
mation height. Here the variations caused by the vari- 
ous PRD approximations are larger. The truncated case 
yields the highest Za, followed by the Doppler case, the 
PRD case, and lowest, the CRD case. As can be seen 
in Table [U the Doppler case approximates the PRD case 
best, with on average a 64 km larger formation height. 
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The differences between the Doppler and PRD cases are 
large where Zg_ is smaU, for example at a; = 4 Mm in the 
lower panel of |6l they are small where is large, such 
as at a; = 1 Mm. 

From the above analysis we conclude that it is possi- 
ble to use Doppler profiles in the Lyman-a and /3 lines 
without significantly infiuencing the formation of Ha. 

5. THE INFLUENCE OF 3D TRANSFER ON Ho 

The conclusions of Sec [2] and Hj namely that it is pos- 
sible to model Ha line formation assuming statistical 
equilibrium and approximating PRD effects in the Ly- 
man lines, makes it numerically feasible to perform 3D 
radiative transfer. Effects from 3D radiative transfer are 
expecte d to be important for st rongly scattering lines like 
Na IDi (jLeenaarts et al.ll2010f l but not so important for 
lines that have a strong coupling of the source function 
to the temperature like the Ca II infrared triplet lines 
(jLeenaarts k, CarlssonI [2009f ). The Ha line is strongly 
scattering as can be seen in the interestin g comparison 
of Ca II 8542 to Ha formation in Fig. 8 of iCauzzi et al.l 
IpoM ). The Ha source function decouples from the 
temperature already in the photo sphere in the ID FAL 
model C of iFontenla et al.1 (|1993[ ) and is afterwards es- 
sentially equal to the angle-averaged radiation field. An 
important corollary of this is that the line-core intensity 
is determined mainly by the photospheric temperature, 
i.e., the line-core resp onse function to temperature peaks 
in the photosphere (jSocas-Navarro fc Uitenbroeld 120041 : 
lUitenbr oek 2006al. 

We therefore solved the radiative transfer problem for 
hydrogen in our entire Bifrost snapshot twice, once with 
full 3D transfer, and once treating each column as a 
plane-parallel atmosphere (hereafter called ID transfer). 
We assumed SE and the approximation to PRD from 
Sec. m No microturbulence was added. The results are 
shown in Fig. [T] 

The top row compares the vertically emergent inten- 
sity between 3D and ID. They show a starkly different 
scene. The 3D computation shows striking, predomi- 
nantly dark, fibril-like features emanating from the pho- 
tospheric field concentrations, which appear to follow the 
magnetic field lines (compare to Fig.[31). They are super- 
imposed on a background of irregularly shaped patches 
of roughly 2 Mm diameter of varying brightness. The 
ID computation instead shows a distinctly photospheric 
scene with granulation and magnetic bright points, much 
more reminiscent of observations in the line wing. 

The second row shows the average formation height. 
The 3D image shows that the structures that appear dark 
in the line-core image are formed highest. The irregular 
background in the intensity image has a lower Za. In the 
ID case the situation is largely the same, except that it 
shows more small scale structure (due to the absence of 
smoothing 3D radiation effects on the opacity) and the 
patches of small away from the magnetic concentra- 
tions. 

The bottom row shows the angle-averaged line-core ra- 
diation field Jj/o at the local optical depth unity. Both 
the ID and 3D images appear very similar to the corre- 
sponding intensity images. This validates the Eddington- 
Barbier (EB) approximation I^oitJ-) ^ Si,o(t = fi) for 
vertical rays {fJ. = 1), and also that the core source func- 
tion is set by scattering (S'j/o(t = 1) ~ Ji^o{t = 1)). 



This also explains why the ID intensity image shows 
granulation and bright points. The Ha source function 
at optical depth unity is given by the Juo{t = 1), which 
in turn is mainly set in the photosphere. In ID there is 
no smoothing of the radiation field, and Juo{t — 1) fii^ a 
given column correlates strongly with temperature at the 
thermalization depth in that column. This temperature 
is lowest in intergranular lanes, higher inside granules 
and highest in magnetic elements. The resulting image 
therefore shows the structure of the photosphere, even 
though optical depth unity lies much higher. 

We conclude that it is essential to include 3D radiative 
transfer effects in the formation of the Ha line core. In 
the next section we analyze the results of our 3D com- 
putation in more detail. 

6. ANALYSIS OF THE 3D COMPUTATION 

6.1. Ha core intensity correlates with average formation 

height 

Fig. [7] suggests a correlation between Ha intensity and 
average formation height for the 3D computation. In 
Fig. IS] we investigate the cause of this effect. The top 
panel shows the probability density function (PDF) of 
the line-core source function at optical depth unity as 
function of emergent Ha intensity in the line core. It 
confirms that the Eddington-Barbier relation is valid. 
The second panel shows the PDF of J^q as function of 
height. There is a rather small spread of J^q around the 
horizontal average of the angle-averaged radiation field 
< Jvo >xy The latter is a decreasing function of height 
up to 2 Mm height, after which it becomes fiat. 

The third panel shows the PDF of the line core 
optical depth unity height as function of emer- 
gent intensity. There is a strong correlation be- 
tween the two quantities between /=4xl0~® and 
/=7xl0~^ergcm~^ s~^ Hz~^ ster~^. At lower in- 
tensity the correlation shows two protrusions with 
a steeper slope. Between /=7xl0~^ and /=8x 
10^^ ergcm^^ s^^ Hz^^ stcr^^ the distribution is bi- 
modal, with maxima at around z = 0.5 Mm and 
z — 1.3 Mm. This bimodality is caused by the mid- 
chromospheric opacity gap: this gap is typically located 
around optical depth unity for those emergent intensities. 
Depending on whether r = 1 is reached above or below 
this gap, a given column in the atmosphere contributes 
to the low or the high cluster of points in the PDF. 

At even higher emergent intensities the r = 1 height 
appears independent of the intensity. The number of 
columns with this intensity is small, so it is unclear 
whether this is significant. The bimodality can be largely 
removed by using the average formation height instead 
of the height of optical depth unity. The former is much 
less sensitive to the location of the t — 1 height than the 
latter. This is shown in the bottom panel. 

The correlation between I and Za (or z(t — 1)) can now 
easily be explained: The Eddington-Barbier approxima- 
tion is valid, and scattering ensures S^o = JvO- The 
radiation field J^o shows relatively little variation with 
horizontal position and its horizontal average is mono- 
tonically decreasing with height. Therefore, larger av- 
erage formation height leads to lower intensity. Above 
z = 1.7 Mm the radiation field is approximately constant 
with height and the formation height and emergent in- 




Figure 7. Comparison of Ha linc-corc formation between full 3D radiative transfer (left-hand column) and treatment of each column in 
the snapshot as a plane-parallel atmosphere (right-hand column). First row: vertically emergent intensity. Second row: average formation 
height. Third row: angle-averaged radiation field at optical depth unity. 
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Figure 9. Comparison of line-core width (top) and temperature 
averaged between t = 0.5 and t = 5 at the wavelength of the 
profile minimum (bottom). 
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Figure 8. Panels, from top to bottom: (1) Probability density 
function (PDF) of the source function at optical depth unity as 
function of emergent Ha intensity in the line-core. The grey line 
indicates S = I. (2) PDF of the angle-averaged radiation field J as 
function of height. The white curve shows the horizontal average 
of J as function of height. (3) PDF of the line-core optical depth 
unity height as function of emergent Hq intensity. The gray curve 
displays the same quantity as the white curve from panel 2, but 
transposed to account for the different arrangement of the axes. (4) 
PDF of the average formation height as function of emergent H« 
intensity. The gray curve is the same as for panel 3. All panels have 
each column individually scaled to maximum contrast for improved 
readability. 




Figure 10. Probability density function (PDF) of the tempera- 
ture averaged between t = 0.5 and r = 5 at the wavelength of the 
profile minimum against the line-core width, after smoothing with 
a 3 X 3 moving boxcar average. The black curve is a quadratic fit 
to the data. 

tensity are no longer correlated. 



6.2. Ha line-core width correlates with temperature 
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Because of the low mass of the hydrogen atom it is 
expected that the width of the Ha Une core is sensitive to 
the temperature in the hne forming region. Because the 
Doppler speed (^2fcT/mH = 13 kms^^ for T = 10 kK) 
is comparable to or larger than typical gas velocities in 
the chromosphere, the core width is not very sensitive 
to velocity gradients along the line-of-sight in the line 
forming region. 

Because we are interested in a diagnostic of the chro- 
mospheric temperature we choose to measure the line 
width close to the line core. For each pixel in our syn- 
thetic image we determined the profile minimum and 
maximum (/min and /max) and then defined the core 
width as the wavelength difference between the intersec- 
tions of the profile and the line / — /min+0.1(/max— /min)- 
We then computed for each pixel the temperature aver- 
aged over the height range between r = 0.5 and t — 5 
at the wavelength of the profile minimum. The height 
range was extended to larger optical depths because one 
looks deeper in the atmosphere at the wavelengths that 
define the line-core width than in the line-core. 

The results are shown in Fig. |9l Comparison of the 
two panels shows that there is a decent correspondence 
between large width and large average temperature, but 
in general no clear correspondence for pixels with low 
width and low temperature. The darker areas in the 
upper panel show a granulation pattern (compare with 
Fig. [3]), while the average temperature panel does not. 
The notable exception are the dark elongated structures 
which clearly show up near the center of both panels. 
Fig.lTOlshows the PDF of the average temperature versus 
the core width after smoothing both with a 3 x 3 moving 
boxcar average. The averaging serves to remove outliers 
in the correlation. 

This figure confirms the impression obtained from the 
images: the two quantities are quite well correlated at 
higher values, but in the low to intermediate regime (be- 
low a line width of 0.6 A) the correlation is rather weak. 

The correlation can be understood using optically thick 
line formation, the Eddington-Barbier approximation and 
the fact that the run of source function with height is 
monotonically decreasing and rather similar for each col- 
umn (cf. Fig. [8]). Because Eddington-Barbier is valid, 
one looks deeper and deeper into the atmosphere as one 
moves away from the line center and one sees a higher 
and higher source function and hence emergent intensity. 
A low temperature means a narrow absorption profile, 
hence steep gradient of the height of optical depth unity 
with wavelength and therefore a narrow line. At higher 
temperature the absorption profile is wider, a less steep 
dz{T = l)/dX gradient and therefore a wider line. The 
actual line width at any given location therefore depends 
on both the temperature, and the gradient of the source 
function with optical depth in the line forming region. 
Variation in the latter is the main cause of the scatter in 
Fig.il 

6.3. No correlation of Ha opacity with temperature in 
the upper chromosphere 

The Ha line opacity is often considered to be sen- 
sitive to the temperature. This is true in the photo- 
sphere and lower chromosphere where the level popula- 
tions are close to LTE, but does not hold in the upper 




T[kK] 

Figure 11. PDF of the relative occupation of the Ha lower 
level as function of the gas temperature at heights between 1 Mm 
and 3 Mm. The black curve indicates the Saha-Boltzmann oc- 
cupation for the average electron density over this height range, 
3.65x10^'' gcm"-^. Each column has been scaled to maximum con- 
trast. 
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Figure 12. PDF of the column mass at the Ha line-core t=1 
height as function of the r=l height. Each column has been scaled 
to maximum contrast. The black curve represents the horizontal 
average of the column mass. 

chromosphere. There the Ha op acity per gram is deter- 
mined by the ionization degree ([Carlsson fc Steinll2002D 
and the radiation field. The ionization degree is insensi- 
tive to temperature variations over time due to the long 
ionization-recombination timescale. The radiation field 
is non-locally determined and therefore independent of 
the local temperature. 

We demonstrate this in Fig.fTTl It shows the n=2 pop- 
ulation relative to the total (neutral plus ionized) hydro- 
gen density for each grid point in our simulation with a 
height between 1 Mm and 3 Mm. At temperatures below 
4 kK the distribution is very wide, at higher temperatures 
the distribution narrows considerably and its peak value 
varies little over the whole temperature range. The black 
curve indicates the LTE value of the relative occupation 
number, which varies by many orders of magnitude. 

The temperature-insensitivity of the upper- 
chromospheric opacity means the opacity is mainly 
determined by the mass density. Therefore, the optical 
depth scale in the upper chromosphere is proportional 
to the column mass. This is shown in Fig. 1121 where we 
show the PDF of the column mass at the t—1 height 
as function of the t=1 height. It shows that if the t—1 
height is above 1.2 Mm, it has a constant column mass 
of roughly 3 X 10~^ g cm^ . 

At lower heights the opacity is sensitive to both the 
temperature and the mass density. The column mass at 
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Figure 13. Image of the line-center intensity. Red arrows indicate 
the direction of the horizontal component of the magnetic field at 
the average formation height. Arrow length is proportional to field 
strength. The contours indicate values of plasma /3 equals 0.01 
(white) and 0.1 (green) at the average f orm ation height. The white 
tick marks indicate the cut used in Fig. 1141 White crosses indicate 
locations of fibrils crossing this cut. This figure is best viewed by 
magnifying the PDF file on a computer screen. 

the lower lying t=1 heights is therefore no longer con- 
stant, and is close to the horizontal average of the column 
mass instead. 

6.4. Fibril-like structures align with the magnetic field 
in a low (3 environment 

Chromospheric fibrils in Ha are associated with the 
presence of photospheric magnetic elements and are con- 
ventionally thought to trace out magnetic field lines, but 
this has so far only been shown for the Ca II 8 542 A line 
(jde la Cruz Rodriguez fc Socas-Navarroll201lD 

Fig. [T3l investigates the alignment of the magnetic field 
in the model atmosphere with the dark fibril-like struc- 
tures in the Ha line-core. The figure shows the emergent 
intensity with the horizontal magnetic field direction and 
the plasma (3 (ratio of gas pressure over magnetic pres- 
sure) at the average formation height overplotted. The 
fibrils are mainly located within the /? — 0.01 contour, 
between the photospheric field concentrations of opposite 
polarity. These fibrils are closely following the magnetic 
field direction. Outside the (5 — 0.1 contour there are 
virtually no fibril-like structures and the image shows a 
more grainy appearance. Note that the semi-transparent 
fibrils at (x, y) = (13.3, 18) Mm (around the upper white 
cross, best seen in the upper-left panel of Fig. [7|) also 
align with the field. 

6.5. Fibril-like structures caused by field-aligned density 

ridges 

We now further investigate the cause of the dark struc- 
tures by showing an yz slice through the atmosphere 
along the cut indicated in Fig. [131 The panels a-d of 
Fig. [T4l show the temperature, mass density, vertical ve- 



locity and magnetic field strength. The fibril locations 
(indicated by the vertical lines) do not show a partic- 
ular correlation with temperature and vertical velocity. 
Panel d again shows they are located in an area with high 
magnetic field and low plasma j3. The mass density panel 
clearly shows that the fibrils are located where there is 
enhanced chromospheric density, which causes the aver- 
age formation height to be larger and thus the emergent 
intensity to be lower (see Sec. 16. ip . 

Panels e-h show the radiative transfer properties of the 
slice. The fibrils appear as local minima in the intensity. 
In the core-forming region the source function is nearly 
equal to the angle-averaged radiation field. 

The contribution functions show a great variation in 
height of formation, exhibit multiple peaks and typical 
extended upward streaks. Away from the regions with 
strong field the average formation height tends to be 
lower, with a strong contribution from around 0.5 Mm 
height. In regions with stronger fields, there is no pro- 
nounced contribution from this depth in the atmosphere. 

Fig. [15] shows the relation between Ha intensity, av- 
erage formation height and chromospheric mass density 
structure. The left panel shows the emergent line-core 
intensity in a subfield of the snapshot that shows fibrils. 
The middle panel shows a 3D rendering of the average 
formation height surface, with the emergent intensity as 
the image on the surface. It shows the dark fibrils as 
ridges with a large average formation height. Finally, 
the right panel shows a 3D representation of the iso- 
surface of the mass density at p = 10~^^gcm~'^, again 
with the emergent intensity as the image on the surface. 
Wherever the average formation height shows a ridge, 
there is a corresponding height increase in the density 
iso-surface. Even though the exact shapes are not iden- 
tical, it is obvious that the dark fibrils between the two 
clusters of opposite polarity field are caused by ridges of 
enhanced chromospheric mass density. The fibrils that 
point away from both magnetic field concentrations are 
also co-spatial with locations of increased mass density, 
but because they are shorter, they appear more as moun- 
tains than as ridges. 

7. COMPARISON WITH OBSERVATIONS 

Fig. [T6| shows observations of a quiet sun region 
taken in the Ha line. The observations were obtained 
with the C Risp Imaging SpectroPolarimeter (CRISP, 
lScharmei~e t al. 2008) at the Swe dish 1-m Solar Tele- 
scope (SST, Schar mer et al.ll2Q03al ) on La Palma on May 
5,2011 at 7:53 UT. The Ha spectral line was sampled 
with 35 line positions, from -2.064 A to +1.290 A, 
and equidistant 86 mA stepping between ±1.290 A. 
The time to complete a full line scan was 8 s. The 
target area was a quiet 60x60"region at disk cen- 
ter centered at {x,y) = (19", 5"). The seeing condi- 
tions were very good and the data q uality benefited 
from t he SST adaptive optics system (jScharmer et al.1 
l2003bD and Multi -Object Multi-Frame B lind Deconvolu- 
tion (MOMFBD. Ivan Noort et all [20051 ) image restora- 
ti on. For more details on the data p rocessing, we refer 
to iRouppe van der Voort et"all (|200l . 

Panel a shows the blue- wing image. It shows granula- 
tion and bright points associated with magnetic concen- 
trations. Panel c shows the line-center image. There 
are two clusters of fibrils centered at (27,7) Mm and 
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Figure 14. Vertical slice through the at mos phere along the cut indicated in Fig. 1131 Vertical black lines indicate the positions of dark 
fibrils marked with white crosses in Fig. 1131 The white or red curve indicates the average formation height, a: gas temperature; b: 
mass density; c: vertical velocity; d: magnetic field strength. The red, green and blue contours indicate plasma /3 equals 1,0.1 and 0.01 
respectively. ; e: emergent line-center intensity f: angle-averaged line-center radiation field; g: contribution function to line-center intensity, 
each column scaled to maximum contrast; h; line-center source function. 
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Figure 15. Display of ridge-like dark fibrils. Left panel: emergent intensity in the central part of the atmosphere. Middle panel: 3D 
surface plot of the average formation height, with the color coding indicating the vertically emergent intensity. Right panel: iso-surface 
of the mass density at p = 10~^^ gcm~'^, with the color coding again indicating the vertically emergent intensity. The dark fibrils follow 
ridges of enhanced chromospheric mass density. 
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(23,27) Mm pointing away from the photospheric bright 
points. There appear to be no fibrils connecting the 
clusters, so it is likely the magnetic field in the im- 
age is largely unipolar. At larger distances from the 
bright points the line-core shows a more grainy appear- 
ance without the elongated structures. This is seen most 
clearly in the upper-left corner and close to the edge of 
the image at the right-hand side. 

The synthetic line-core image in the upper left panel 
of Fig. [7] shows a similar structure, it shows fibril-likc 
structure close to the photospheric magnetic concentra- 
tions and grainy structure elsewhere. The simulations 
lack the dense packing of the fibrils in the observations 
and shows typically shorter fibril-like structures. The 
width of the dark fibrils in the observations are similar 
to the width of the dark fibril-like structures in the sim- 
ulations. 

Panel b shows the line-core width, determined in the 
same way as for the simulations (see Sec. l6.2p . The fibrils 
show largest width. The areas with lower width do not 
show a granulation pattern, instead it shows a pattern 
reminiscent of reversed granulation or acoustic shocks. 
Comparison of panel b and Fig. El) shows that the syn- 
thetic fibril-like structures show up much less in the line- 
width measurement than the observed fibrils. The signa- 
ture of the granulation pattern in the simulations shows 
that the simulated internetwork has a lower Ha opac- 
ity than the sun. The range of line-core widths in the 
simulation is very similar to the one measured in the ob- 
servations. 

Panel d finally shows the gas temperature derived from 
the core width assuming the fit from Fig. (TU] It shows 
temperatures between 4 and 8kK for internetwork areas. 
The fibrils have a higher temperature of around 10 kK. 
Note that the brightest patches in the line-core image 
do not correspond to the largest core width or derived 
temperature. 

8. DISCUSSION AND CONCLUSIONS 

We have presented results of a radiative transfer com- 
putation of the Ha line based on a model atmosphere 
computed with a state-of-the-art radiation-MHD code. 
This is the first time that a simulation reproduces chro- 
mospheric fibril-like structures in synthetic observables. 

The crucial ingredient for the formation of fibril-like 
structures is the full 3D evaluation of the transfer equa- 
tion. The Ha source function decouples from the temper- 
ature already in the photosphere and the source function 
at optical depth unity in the chromosphere is therefore 
set by radiation emitted in a large volume and hence 
fairly smooth. 

In our simulation the emergent line-core intensity cor- 
relates with the average formation height: the lower the 
intensity, the higher in the atmosphere the photons are 
originating. It is a very different explanation than the 
one that is typically assumed by cloud modeling, where 
variation in emergent intensity is assumed to be caused 
by varying thermodynamic properties in a cloud that 
modulate the intensity emitted by the photosphere be- 
low. 

We expect that this intensity-formation-height corre- 
lation is also valid on the sun, as it is based only on the 
assumption that the Ha line is strongly scattering within 
its formation height range, so that the source function is 



set by scattering and the angle-averaged radiation field 
shows relatively little variation with horizontal spatial 
position. However, the shape of the correlation curve de- 
pends on the exact run of the source function with height. 
The latter depends on the temperature at the thermaliza- 
tion depth and the chromospheric opacity. These will be 
different in, for example, internetwork, network or plage 
and different again in our simulation. When comparing 
intensities for small areas on the sun one can conclude 
that the lower the intensity the higher in the atmosphere 
one looks, but intensities in very different regions cannot 
be directly compared. 

A largely similar argument holds for the correlation 
between line-core width and gas temperature. This cor- 
relation can be simply explained by optically thick line 
formation and the Eddington-Barbier approximation and 
is therefore doubtlessly valid on the sun too. The ex- 
act shape of the correlation depends on the variation of 
source function and opacity with height. The fit shown 
in Fig. [To] is therefore likely not universally valid. 

In addition, the method is sensitive to velocity gradi- 
ents that are of the same order as the line width. Such a 
gradient causes an increase in line width compared to 
an identical atmosphere without a velocity field. For 
example , the asymmetric profiles of 'rapi d blueshifted 
events' (jRouppe van der Voort et all I2009D would show 
up as having increased temperature using the core- width 
measure, without this necessarily being the case. 

Still, if we assume it is valid and apply it, our 
observations yield internetwork temperatures compat- 
ible with radiation- MHD simulations of the internet- 
work chromosphere (ICarlssori fc SteinI [l995t iSkartlienI 
120001: iWedemever et al.l 120041: iLeenaarts et a l. 2007). 
We find temperatures of 10 kK for Ha fibrils. This 
is comparable to values found in many studies that 
use cloud modeling (e.g., 'Heinzcl & Schmicdcr 1993; 
iTsiropoula fc SchmiedeH[l997,: Bostanci 2011, ) However, 
all cloud modeling requires addition of a 'microturbulent 
velocity' of typically 5 to 15kms~^ as an ad hoc param- 
eter to derive the temperature from the observed cloud 
model parameters. The profiles in our simulation have 
similar widths as the observations and yield similar tem- 
peratures without the need of invoking microturbulence. 

The line formation properties of Ha show a dou- 
ble character. At heights below 1 Mm its opacity is 
temperature-sensitive because of its excited lower level. 
At larger heights this temperature sensitivity disappears. 
There the opacity is instead mainly sensitive to den- 
sity, much more akin to resonance lines, which have the 
ground state as a lower level. 

Our simulation supports the commonly held notion 
that fibrils indeed trace the magnetic file. More precisely 
we show that the orientation of the fibril-like structures 
in the synthetic Ha core image align with the horizontal 
component of the magnetic field at the average formation 
height computed for each pixel separately. This supports 
the practice of using Ha imagery to constrain and test 
the qua lity of field extrapolation s (e.g ., I Judge et al.ll2010l : 
iJing et al. 2011; Rea rdon et a l. 20111 ). 

The dark fibrils-like structures in the simulation are 
caused by ridge-like density enhancements in the chro- 
mosphere that lead to higher formation height and hence 
lower emergent intensity. It is unlikely the density ridges 
are the only way to create fibrils on the Sun. The sim- 
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ulation has a limited resolution and spatial extent and 
still lacks treatment of physical processes that might be 
important in the chrom osphere, such as hea ting due to 
neutral- ion drag (jKrasn osclskik h et alJl20l"ol ). 

We can now answer the question what makes Ha such 
a good line to observe the magnetic structure of the chro- 
mosphere. The main reason is that it happens to have so 
much opacity in the line core that it almost always forms 
in the low plasma beta regime where the magnetic field is 
the main structuring agent in the atmosphere. The low 
mass of the hydrogen atom makes the line wide, so that 
the fixed-wavelength line-core intensity is only weakly 
modulated by the velocity field. Instead, the correla- 
tion of formation height and intensity and the constant 
column mass at the t=1 height make the line-core in- 
tensity a tracer of the chromospheric mass density. It is 
the variation of the mass density caused by the magnetic 
field, waves and shocks that gives rise to the dramatic 
structures seen in Ha. 
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